Relaxation, Emission Modes and Limit 
Temperatures in Ni+Ni HIC. 

Armando Barranon *; J. A. Lopez C. Dorso. * 
November 15 th 2003 



Abstract 

A dynamical stability analysis is performed for Ni+Ni central collision 
at intermediate energies, showing that in chemical, thermal and dynamical 
equilibriums are reached at an early stage of system evolution. This is ob- 
tained by computing the relaxation times of the quadrupole momentum, 
speed of sound and electric charge density of the system. This way, frag- 
ment emission modes at low and high energies as well as the qualitative 
behavior of the limit temperatures are determined. 

1 Introduction 

In a heavy ions collision, a compound nuclear system is initially formed 
that will expand as a consequence of the thermal or compressive pressure. 
Expansion leads the nuclear system into the instability region where the 
speed of sound attains an imaginary value. At this instability region, sys- 
tem is disassembled into fragments and nucleons. Fragmentation results 
from an adiabatic expansion of hot nuclear matter produced in the col- 
lision that leads the system into the mechanical instability region, where 
the system fragments into clusters and nucleons due to the increase of 
density fluctuations 0. Neutron rich Heavy ions collisions lead to tran- 
sient states of nuclear matter with high spin asymmetry and with a great 
thermal and compressive excitation, allowing the exploration of the prop- 
erties of nuclear matter whose neutron contain goes from pure neutron 
nuclear matter up to symmetric nuclear matter Mller et. al. has 
shown that the liquid-gas phase transition in asymmetric nuclear matter 
is of second order instead of the first order phase transition observed in 
symmetric nuclear matter [3]. 

Several important aspects of nuclear collisions depend on the isospin, 
such as the nuclear stopping, the collective nuclear flow as well as the bal- 
ance energy and pre-equilibrium nucleon emission. Bao An-Li et. al. have 
found that the mean nuclear field domains completely the density fluctu- 
ations as well as the isospin fluctuations, meanwhile two-body scattering 
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meaningfully influences in the subsequent isospin growth. Both fluctua- 
tions are increased in the dynamically instable systems meanwhile in the 
chemically instable systems only the density fluctuations are increased. 
The magnitude of both fluctuations diminishes when increases the isospin 
asymmetry resulting from the decrease of the isoscalar attractive mean 
field, due to the increase of the neutron symmetric repulsive potential in 
neutron rich nuclear matter [I]. Beaulieu et. al. analyzed experimental 
results of heavy ions collisions and found long relaxation times for surface 
fragment emissions and short relaxation times when fragments are emit- 
ted from the bulk 0. In the soft region of the state equation, the pressure 
change is much smaller and hence a small speed of sound is attained since 
it depends directly on the pressure energy gradient: V 2 = ^ . 

In this soft region, the collective dynamics of dense and hot matter 
formed in heavy ion collisions can be determinant. As a matter of fact, 
a small speed of sound delays the expansion of compressed matter and 
leads a lighter transverse collective flow [Sj. As shown via hydrodynamic 
models, heavy ion collisions reach the soft region of the state equation 
when the incident energy varies, building up a fireball that lasts a long 
time, leading to a minimum in the energy dependence of the collective 
transverse flow 0. 



LATINO model [8] uses semi-classical approximation to simulate inter- 
mediate energy heavy ions collisions and replicates binary interaction via 
Pandharipande potential. This potential is built up with a linear combi- 
nation of Yukawa potentials, whose coefficients are fitted in order repro- 
duce the properties of nuclear matter. Clusters are recognized using an 
Early Cluster Recognition Algorithm, that optimizes the configuration in 
the energy space. Ground states are produced generating a random con- 
figuration in phase space, gradually reducing the speed of the particles 
confined into a parabolic potential, until the theoretical binding energy is 
reached. Projectile kinetic energy varies in a range going from 800 MeV 
up to 2000 MeV in the center of mass system, with 200 replicas for each 
energy. Numerical integration of the equations of motion is performed 
with a Verlet algorithm with time intervals that ensure energy conser- 
vation up to 0.05converge to a power law in this range of intermediate 
energies 



In the Minimum Spanning Tree in Energy Algorithm (MSTE) , a given 
set of particles i,j,...,k belongs to the same cluster d whenever: Vi £ d, 
there exists j £ Cj/ey < where ey = V(ry) + (pi — p-,) 2 /4/x , and fi is 
the reduced mass of the pair i,j. 

Fragment formation time, tff, can be defined as the time when the 
system breaks up in a defined way, namely after tff fragments evaporate 
only some monomers. In order to estimate tff the Microscopic Persistence 
Coefficient P is used, |10| : 
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Latino Model. 



3 Methodology 




(i) 
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where X = {C;} and Y = {C|} are two partitions, bi is the number of 
particle pairs belonging to cluster d of partition X, a; is the number of 
particle pairs belonging to cluster d and also belonging to a given cluster 
C'j of partition Y, rii is the number of particles in cluster d . 
It is convenient to study the time evolution of the quantities: 

P + [X(t)] = (P [X(t),X(t - oo)]) co „ isiona (2) 
P- [X(t)] = (P [X(t - co), X(t)])„ niaima (3) 

P dt [X(t)] = {P[X(t),X(t + dt)}) comsions (4) 

where X(t) is a partition at time t, X(t — > oo) is an asymptotic partition, 
and represents an average on all the collisions, is a partition identical to 
the partition , except for the fact that a nucleon has been evaporated in 
each cluster. 

x = {n/p)\ y>0 / ( n /p)\ y <o ( 5 ) 
where (n/p)\ >0 and (n/p)\ <Q are the ratios between the number of 
neutrons and protons moving forward and backwards, respectively. Dy- 
namical instability is set when the squared speed of sound is negative. 
Adiabatic speed of sound is given by: 

VS = (l/m) [(dP/dp))s] (6) 

VS = (l/m) [(10/9) < E K > +a(p/po) + b{p/p Q ) a ] (7) 

where < Ek > is the mean kinetic energy per nucleon, a=-358.1 MeV, 
b=304.8 MeV and a — 7/6 are the parameters of the soft state equation. 
When Vq < nuclear matter is instable with respect to density fluctu- 
ations, leading to dynamical instabilities Heavy residuals thermal 
equilibrium can be examined studying the quadrupolar momentum Qzz 
defined by: 

Qzz = J drdp(2nf [2pt - p 2 y p 2 x ] f(t, 7, t) (8) 

where f{r,~p*,t) is Wigner function. Clearly Qzz — is a necessary 
but not sufficient condition for thermal equilibrium 1121 . Mean Velocity 
Transfer is given by: 

MVT, = J2 \Vu(t + dt) - Vki(t)\ (9) 

4 Relaxation and Instabilities. 

Collision Disordered Mode is the only one present in the uniformly excited 
systems, and is responsible of fragment production as well as the exterior 
flow dispersing the fragments 1131 1141 . MSTE can be used to study biggest 
fragment size, total multiplicity and persistent coefficients. 

Figs, la, b and c show time evolution of these three quantities. Due 
to the initial violent collision, some particles attain lots of energy and are 
unable to be bound to a fragment even when these particles might be 
inside or close to a fragment. This reduces the size of the biggest MSTE 
cluster in this early stage and increases the total multiplicity. Fig, 2 shows 
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Ni+Ni collision stages for an energy of 1400 MeV. The first stage ends 
up with a peak of kinetic energy transported by the promptly emitted 
particles, KE(PEP), signed by the peak of the wide continuous line. The 
second stage finishes up with the attenuation of the intermediate fragment 
production (IMF) , shown as a peak in the continuous curve. Between 
these two peaks, another peak is seen in the intermittent wide curve of 
the Mean Velocity Transverse (MVT), signing the onset of the instabili- 
ties. As can be seen, MVT produce instabilities leading to light particle 
emission and subsequent intermediate fragment emission (IMF). Once the 
continuous curve is attenuated, a final stage is installed with slow light 
fragment emission and system freezing-up. Persistent coefficients P + and 
P~ computed using MSTE partitions, can be sed to understand how the 
partition attains its microscopic composition at an asymptotic time. 

In this study dynamical, thermal and chemical instability signals are 
used to detect the relaxation times for several projectile energies, that are 
expected to reach a peak in the second collision stage and attenuate at the 
third collision stage. Fig. 3 shows the time evolution of the parameter 
A with a peak signing the entrance into the chemical instability zone, 
for the Ni+Ni central collision with a projectile energy of 1400 MeV. 
When the system attains chemical equilibrium, the parameter Areaches a 
value close to zero. Fig. 4 shows the time evolution of the normalized 
quadrupole momentum with a peak signing the entrance into the thermal 
instability zone, for the central collision Ni+Ni, for a projectile energy of 
1400 MeV. When the system attains thermal equilibrium, the normalized 
quadrupolar momentum reaches a value close to zero. Fig. 5 shows time 
evolution of the squared adiabatic speed of sound with a peak signing the 
entrance into the dynamical instability zone, for Ni+Ni central collision 
with a projectile energy of 1400 MeV. When the system reaches thermal 
equilibrium, the squared adiabatic sped of sound reaches a value close to 
zero. 

5 Conclusions. 

Relaxation times obtained via these three signatures show that when the 
projectile energy is slow, the most important relaxation time is the ther- 
mal one, as expected since the available energy is spent in deformation 
and fragments are emitted from the surface (Fig. 6). But when the pro- 
jectile energy is high, fragments are emitted from the bulk, and the three 
relaxation times are the same. Therefore, we expect to obtain low tem- 
peratures when the residual is big, since typically this case corresponds 
to low projectile energies and low excitations. Fig. 7 shows the time evo- 
lution for a typical experiment of the central collision Ni+Ni with a low 
projectile energy equal to 1100 MeV, where as can be seen, the compound 
suffers deformations up to an asymptotic state and a long quadrupolar 
relaxation time is expected. In a similar fashion, we expect high tem- 
peratures for a high excitation and a high projectile energy, since the 
relaxation times are short and the system disassembles almost immedi- 
ately, spending the available energy on heating the system and leading to 
higher temperatures. When the projectile energy is increased, fragments 
are no longer emitted from the surface but from the bulk emission, and 
the dynamical behavior of the system is expected to change, providing dif- 
ferent limit temperatures and distinct residual sizes. Work supported by 
the National Science Foundation (PHY-96-00038), Universidad de Buenos 
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Aires (EX-070, Grant No. TW98, CONICET Grant No. PIP 4436/96), 
and Universidad Autonoma Metropolitana-Azcapotzalco (Laboratorio de 
Computo Intensivo). 
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